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The dynamics of a kicked quantum mechanical wavepacket at a quantum resonance is studied 
in the framework of Floquet analysis. It is seen how a directed current can be created out of a 
homogeneous initial state at certain resonances in an asymmetric potential. The almost periodic 
parameter dependence of the current is found to be connected with level crossings in the Floquet 
spectrum. 

I. INTRODUCTION 

The kicked rotor is, due to its relative simplicity, one of the most studied and best understood models of chaotic 
mechanics. In essence, it models a particle that is exposed to a sinusoidal potential during periodically repeated, 
delta-function shaped kicks. The quantum mechanical version of the kicked rotor has attracted special attention 
lately, due to the realization of the model in optical lattices lii, i2j l3| • It turns out that there exist two peculiar 
quantum mechanical effects that distinguish the quantum kicked rotor from its classical counterpart. These two, 
mutually exclusive, effects are Anderson localization in momentum space, which occurs for generic (almost all) values 
of the kickingperiod 0], and quantum resonances, which occur when the kicking period matches a resonance criterion 
0, 0, IS IE I3;H; 13 ■ The latter phenomenon, quantum resonances, is the subject of this paper. 

Classically, the energy growth of the kicked rotor is diffusive, i. e. linear. In the quantum case and off resonance, 
the energy growth will also be diffusive but eventually saturate due to Anderson localization in momentum space. 
In contrast, a quantum resonance will result in an energy growth that is quadratic in time. Such resonances occur 
when the period of the potential kicks on the particle matches a rational value times the so-called Talbot time of the 
system. In the context of optical lattices, the corresponding frequency is identical to the recoil frequency 10]. It is 
common to state the resonance criterion in terms of an effective Planck constant which is defined as a combination of 
the physical parameters. 

We have recently found that a slight generalization of the kicked rotor to an asymmetric, sawtooth-shaped potential 
will, combined with quantum resonances, result in a ratchet effect of sorts - a directed current which increases linearly 
with time p^ . In general, a classical or quantum particle in a flashing periodic potential will not pick up a finite 
velocity even if the potential is asymmetric, i. e., there is no ratchet effect. However, quantum resonances may change 
the situation and allow for a directed current. The effect hinges on the fact that a specific momentum eigenstate is 
chosen as the initial state, since a proper averaging over all of phase space would necessarily cancel out any directed 
current. Nevertheless, it was seen in Ref. that the proposed setup could give a constant acceleration to an initially 
zero-momentum plane wave. 

The spectral properties of the quantum kicked rotor were first studied by Izrailev and Shepelyanskii [H 1^1- It 
was found that the quasienergy spectrum is in general discrete, but in the case of a quantum resonance it forms a 
band structure with a continuum of quasienergies. This was found to explain the quadratic energy growth. In the 
present paper, the same type of analysis is exploited in order to understand how a directed current is created out of 
a motionless initial state with the help of an asymmetric potential. In Sec. ^we describe the problem and set up the 
definitions. In Sec. Illll we show how the general solution of the problem at hand is constructed and cover the cases of 
two simple resonances. In Sec. II VI a resonance of special interest is analyzed by perturbative and numerical means. 
Finally, in Sec. we summarize and conclude. 
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II. PRELIMINARIES 



We wish to study the Schrodinger equation 

n— 1 

This equation describes a particle subject to a periodic potential Uix) that is flashed on for short, periodically 
repeated pulses. It can be realized with atoms in an optical lattice 0, Si 01 ; the kicked rotor problem corresponds 
to the case U{x) oc since, while the superposition of two sines has been predicted to result in a directed current [lol |. 
The parameter S is a dimensionless, effective Planck constant. Units are chosen so that the spatial periodicity of the 
potential is 27r and the temporal period of the flashing is unity. 
The time development of the wave packet is given by 

/•2t °° 2 

ij{x,t + l)=U[tP{x,t)]^e-'^'^''^ I dx' e~'^^''-''">-'^^i^{x',t), (2) 

k— — QO 







where the operator U was implicitly defined. We assume that the wave function has the same spatial periodicity as 
the potential, and hence the momentum variable k is restricted to integers. The extension to non-integer momenta 
is straightforward but not important to the objectives of this paper. For notational convenience we have defined 
V{x) = Uix)/k. 

The time development can be written in terms of the Floquet states Wj (x) , which are eigenstates with associated 
quasi-energies loj of the evolution operator for one temporal period: 

e-^-^w,{x)=U[w,{x)]. (3) 

Since U is unitary, the quasi-energies ujj are real. The time evolution of an initial state "0(2;, 0) is constructed as 
follows: 

i^{x,t) = Y,c,e-^'^^'wj{x), (4) 



with 



dxWj{x)*^|J{x,0). (5) 



Note that since the evolution operator U propagates the system for one unit of time At = 1, Eq. is valid for integer 
t only. 

This type of analysis was first performed for the kicked rotor in Refs. [^0. It was found that quantum resonances 
are associated with a banded quasienergy spectrum. In the present paper, the same type of analysis will be employed 
with a specific goal in mind: It was found numerically in Ref. 10] that when the effective Planck constant k takes on 
values that are integer or half- integer multiples of tt, the ensuing resonant behavior may result in a directed current. 
We shall now investigate how this is reflected in the quasienergy spectrum. 



III. FLOQUET STATES AT RESONANCES 

This section contains mainly a brief repetition and slight generalization of the findings of Ref. 0], so that the 
formalism can be applied to the case of an asymmetric flashing potential. In order to get acquainted with the system 
and define the terminology, we start with the simplest case, where the Planck constant is an integer multiple of An. 
It is well known that in this case, the density stays unchanged at all times, the mean momentum increase is zero and 
the energy increases quadratically with time for any potential V{x) ^M,M^- 

When k — Att, then exp(— i^/c^/2) = 1 for all k, and the unitary operator U in Eq. jSJ simplifies to 

k— — OQ 
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Already from here we can see the solution to the full problem, but in order to prepare ourselves for more complicated 
cases we solve for the time evolution using Floquet analysis. The eigenvalue equation reads 

e-'^^''^Wj{x)^e~'''^w,{x). (7) 

The phase factor on the left-hand side is space dependent, but that on the right-hand side is not. This has a solution 
if Wj{x) is nonzero only for a discrete number of spatial points x. The discrete index j has to be changed into a 
continuous one, xq, and the solution for the continuous set of eigenstates, {wxo}o<XQ<2Tr, is found to be 

Wxoix) ^ 5{x - xo)j i^xo =y[xo). (8) 

The expansion coefficients for the initial state are 

dxWxo{x)*ij{x,Q) = tp{xo,Q), (9) 

and the wave function at integer time instances t is 

^(x,t)= /da;oc,„e-*"-o^^^(a;) =^(a;,0)e-^*^(^). (10) 



We conclude that the system remains unchanged at all times except for a multiplicative phase factor; there is no 
transport. 

We are now prepared to discuss the more general case h = 7r/m, where m is an integer. (The case ^; = 27r is more 
easily solved by simpler means |10].) We have for the sum over momenta 

oo 2m — 1 

E e-"^e-'=-= (11) 



where 



2m-l 

I 

A 

2m 

n=0 



2m-l 

_ e-*^("'"2"^). (12) 



The Floquet equation becomes 



2m-l 

„(a:)=e-^(-) J] A,Wx,,{x - j-). (13) 

^ — ' TO 



J=0 

By the same argument as above in the k = Att case, the solution must be of the form 

2m- 1 



oi^) = y2 cnxoiSix- xo- I — )- (14) 



TO 

1=0 



We have come to the important conclusion that at a quantum resonance, each Floquet eigenstate is nonzero only 
at a discrete set of points. As a result, the Floquet spectrum is divided into a number of discrete bands, within 
each of which the quasienergy depends on the continuous parameter xq ■ Now, inserting this ansatz and equating the 
coefficients yields the equations for the factors a^o,; 

2m-l 

e"'""oaxo,n= E ^^"'"^0,^ (15) 

where 

Mr^i = e-'^("°+"-/™U„_z. (16) 

The indices I and n are modulo 2to. There are 2to solutions to this eigenvalue equation, which will be labeled 
a^gi, with /i = 0..2to. But there is a degeneracy, since we can choose phases such that w'^^-^^rTr/m^^) ~ '^xoi^)' '^^ 
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xq+ttt /rn,l ^Xo,l+r- 



This degeneracy needs to be taken care of by restricting < a;o < 7r/m, in order to avoid 
overcounting. 

Let us now use these eigenstates to construct the long-time development of the wavepacket. We specialize to the 
case of homogeneous initial conditions, '^{x,{)) = l/\/27r, whereby the coefficients c become 



2m-l 



(17) 



and 



(18) 



where Ix is the integer that fulfills nlx/m < x < ^{Ix + l)/m. After the completeness of the basis w^^ has been 
exploited, we are free to extend the definition of the amplitudes a to the whole range < xq < 27r, and there results 



(19) 



The goal is to calculate the time dependence of the momentum. It is given by 



dx 



(20) 



The sum can divided into a part that oscillates in time and a constant part. The constant part will be responsible 
for a linear increase of the momentum, if there is one. This term is denoted by F and is given by 



1 /■27r 

1,1' 



(21) 



This is the final expression for the rate of momentum increase for an initially homogeneous wavepacket. 

The solution method outlined above will now be applied to the case S = tt, i. e. m = 1. We have exp(— i7rA;^/2) = 1 
for even k and — i for odd fc, and the kernel is 



E 



^-ikx^-M- = ^6{x) + ^S{x - tt). 



(22) 



The solutions to the Floquet equation consist of two bands, whose discrete indices ^ will be labeled + and -. The 
eigenvectors are 



-iAV 



It 



sin AT/ 



2 V' Vl + sin^AF, 
^ / ^ sinAF 



Vl + sin^ AV, 



1/2 



1/2 



(23) 



with the eigenvalues 



1-i 



cos AF ± i Vl + sin^ AV 



(24) 



where AV = {V{xo) - V{xo + 7r)/2, and V = {V{xo) + V{xo + tt)/2. Thus, 



uJxo = -J + V{xo) =F arctan 



1 + sin^ AV{xo) 
cos Ay(a;o) 



(25) 
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Now insert this into the general expression, Eq. (|21|l . for the time development. There results for the time- 
independent part of the force 



dx 



dV{x) 
dx 



-(IT 



piAV „ 

V(i± 



It 



dxV'ix) 1 - 



1 + 



(26) 



where s and c are short for sin Ay (x) and cos AV^(x), respectively. The presence of a term odd in AV{x) is crucial. 
Now use the periodicity of V to transform 



dxV'ix) 



1 



-V'{x)-^ + V'{x + tt)- 



27r 



2ir 



1 + S2 
SC 



dxV'ix) - / (ia;Ay'(a;)da;- 

' 1 + 

sin AV^ cos AF 

a;=o Jx=o 1 + sin^Ay 



dV{x) 



d{AV) 



= 0, 



(27) 



by the periodicity of the potential. This could be done because the integrand of the second term depends on the 
coordinate solely through the potential difference AV^. This concludes the demonstration that the drift is zero in an 
initially homogeneous system for any potential V at the resonance k — n. 



IV. RESONANCE AT = 7r/2 

The half-integer resonances are especially interesting since they are known to result in directed transport even if 
the initial state is homogeneous 0|. We therefore study the case k — it/2 with special care. The momentum sum is 
now 



^ 2 



e-'^'^Six) + 5{x - 7r/2) - e-'^^'^Six - tt) + 5{x - 37r/2) 



(28) 



and correspondingly the matrix M for the Floquet eigenvectors reads 



M = V- 



1 



/ e-'-''^ 1 -e 



-i7r/4 



-i7r/4 -j^ \ 

1 _e-»'^/4 



\ 1 



g-j7r/4 ^ g-i7r/4 



(29) 



where 



(30) 



We first solve this system by perturbative means. The quasienergies and eigenvectors to zeroth order in the potential 
V{x) are 







= 71", 


"(0) = ( 






= 0, 


"(0) = 






7r/4, 


"(0) = 


^(0) 




7r/4, 


«(0) = 



(31) 



so the correction to first order is 



«(i) = -i{-Vo + Vi-V2 + V:,){lA,lAf + 



-an/8 



V2T72 



(Ao,-Ai,-Ao,Ai)^ 



6 



u- -2 




FIG. 1: rate of momentum increase on an initially homogeneous wavefunction for a = 0.3, as a function of potential strength 
U, assuming the form of Eq. 13311 for the potential. 



1 g»'r/8 
-l{-Vo + V1-V2 + V3){1, -1, 1, -if + ^ 

2 V2W2 



(-Ao,-Ai,Ao,Aif , 



4) 



-\/2Ai(f 



i7r/4 



l,e*"/^l) 



1) = -^/2Ao(l,e'-/^l,e^-/T, 



*(i) 



(32) 



where we defined Vj = V{x + j7r/2) for j — 0, 1, 2, 3; Aq = Vq — V2; and Ai = Vi — V3. These eigenvectors can now 
be inserted into the expression for the force, Eq. (|21|) . In order to obtain a closed expression, one has to make an 
assumption about the potential. We make the physically motivated choice 



V{x) = C/(sinx + asin2x), (33) 

which models an atom subject to standing laser waves flfll- There results for small U 

F ^ {3 + 2V2)aU^ - {l + V2)a^U^. (34) 

In order to go beyond perturbation theory, the k — tt/2 problem has to be solved numerically Again we assume 
the form Eq. (|33|) for the potential. Figure ^ depicts the numerically calculated force F as a function of potential 
strength U, for the choice a = 0.3. The force rises initially as the third power of U as the perturbative calculation 
indicated, but for longer times the time dependence displays an approximate periodicity of F with period tt as a 
function of U. We shall now see that it can be related to the structure of the Floquet spectrum. In Fig. [3 we display 
the numerically obtained eigenvalues uii^^ for a range of values of U. It is seen that a level crossing as a function of the 
continuous parameter xq appears when U = 3.3 and then separates when U is increased. At this point, two energy 
bands momentarily merge into one. This crossing is reflected in Fig. ^ as a pointed feature in the force curve F[U). 
Closer inspection reveals that there is a discontinuity in the first derivative at the turning point, but there does not 
appear to be a cusp. 

Level crossings in Floquet s pec tra have been seen in various contexts to be associated with resonance phenomena 
[TTIIT^ . However, in Refs. [Tllll2j. the situation is different: the eigenvalue spectra are discrete and the level crossings 
occur as a control parameter is varied. In the present system, the eigenvalues form a continuous set and the level 
crossings take place in the two-dimensional space formed by the continuous index xq and the parameter U; in other 
words, the crossing is a merging of two quasienergy bands. The crossings are not easily interpreted as a resonance 
phenomenon, but they are connected with the sharp turning points of the force as a function of U that can be 
seen in Fig.^ In fact, every such sharp turning point coincides with a level crossing, i. e., a merging of bands, in the 
continuous spectrum. In Fig. 3, we see how a number of level crossings appear around U — 19. Here, the minimum 
in the curve F(U) is smoother, and accordingly it is associated with a sequence of level crossings instead of just one. 

A straight level crossing originates from two eigenvectors that have different symmetry in the sense that the 



FIG. 2: Floquet eigenvalues tj^^^ as functions of xq for four choices of potential strength U. There appears a level crossing at 
X fa 0.4 for U — 3.3 which disappears again for larger U; this signals a sharp feature in the curve for the force F in Fig. 
Because of periodicity, only the range 0<a::<7r/2is displayed. 




0.1 0.2 0.3 0.4 0.1 0.2 0.3 0.4 



X X 

FIG. 3: Floquet eigenvalues cul^g as functions of xq for four choices of potential strength U. 

overlap with respect to a small change 6U of the control parameter U vanishes; 

ia^'f—SUd^'=0. (35) 

(Alternatively, a small variation in the index a;o can be considered, since the level crossing takes place in the two- 
dimensional space of U,xo.) The generic situation is that the overlap is nonzero, which results in an avoided crossing 
and a mixing of the eigenvectors; when the overlap vanishes the energies cross and the eigenvectors do not mix. 
Indeed, this decoupling of eigenvectors is clearly seen in the plots of the components a of the eigenvectors, a selection 
of which is displayed in Fig.^ The avoided level crossing is accompanied by the crossing at a; « 0.9 of the components 
of the corresponding eigenvectors. Close to the turning point the slope gets steeper and eventually the two lines do 
not cross at all, indicating that the overlap between the modes vanishes. This is, in turn, reflected in the force integral 
(|21|l and manifests itself in the sharp bending of the curve. 
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U=3.1 U=3.2 




0.5 0.6 0.7 0.8 0.9 1 0.5 0.6 0.7 0.8 0.9 1 



U=3.3 U=3.4 




0.5 0.6 0.7 0.8 0.9 1 0.5 0.6 0.7 0.8 0.9 1 



FIG. 4: A selection of eigenvector components qJ^j as functions of x, for a few choices of U. 

V. CONCLUSIONS 

We have investigated the Floquet spectrum for a wavepacket subject to periodic forcing, with special attention to 
quantum resonances, where the effective Planck constant is equal to a rational multiple of tt and the quasienergies 
form a band structure. We derived an expression for the force on the wavepacket in terms of Floquet eigenstates and 
concluded that a nonzero mean velocity is obtained from homogeneous initial conditions only at minor resonances 
when the Planck constant is equal to a half-integer multiple of tt. The oscillatory dependence of the current on the 
potential strength was investigated for the special case of a potential composed of two harmonics, and it was seen 
how turning points in the curve of drift as a function of potential strength arise from level crossings in the Floquet 
spectrum. 
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